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Abstract 

We consider the classical and quantum dynamics in M(atrix) theory. Using a 
simple ansatz we show that a classical trajectory exhibits a chaotic motion. We 
argue that the holographic feature of M(atrix) theory is related with the repulsive 
feature of energy eigenvalues in quantum chaotic system. Chaotic dynamics in 
N = 2 supersymmetric Yang-Mills theory is also discussed. We demonstrate that 
after the separation of " slow" and " fast" modes there is a singular contribution from 
the "slow" modes to the Hamiltonian of the "fast" modes. 
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Recent advances in string theory have led to the discovery of dualities between five 
known superstring theories. These theories are expected to be obtained by taking vari- 
ous limits of the conjectured eleven-dimensional M-theory [0, |^. At low energies/large 
distances M-theory is described by eleven- dimensional supergravity. Banks, Fischler, 
Shenker and Susskind have proposed that M-theory in the infinite momentum frame 
is described in terms of a supersymmetric matrix model, the so called M(artix) theory. 
Moreover, the only dynamical degrees of freedom are Dirichlet zero-branes and the calcu- 
lation of any physical quantity of M-theory can be reduced to a calculation in the matrix 
quantum mechanics. A system of N Dirichlet zero-branes is described in terms of nine 
N X N Hermitian matrices Xi,i = 1,...,9 together with their fermionic superpartners. 
The action can be regarded as ten-dimensional SU{N) supersymmetric Yang-Mills theory 
reduced to (0-1-1) space-time dimensions: 

S = J dtTr{^DtX,DtX, + ^[X,, Xj][Xi, X,]) + (ferrmons), (1) 

where = dt + iAq. The action (|l]) was considered in the theory of eleven-dimensional 
supermembranes in |^, ^ and in the dynamics of D-particles in J^, 0. In the original 
formulation of the conjectured correspondence between M-theory and M(atrix) theory 
the large N limit was assumed. A more recent formulation [^] is valid for finite N. The 
M(atrix) theory is interesting because it could provide a non-perturbative approach to 
quantum gravity. Therefore it is important to investigate exact properties of the model 

In this paper we study the bosonic part of the dynamical system (|l|). We show a 
complicated chaotic behavior of classical trajectories and discuss its quantization. The 
appearance of chaos in a classical system means that we cannot trust to the ordinary 
perturbative analysis of the corresponding quantum system. We demonstrate that after 
the separation of "slow" and "fast" modes there is a singular contribution from the "slow" 
modes to the Hamiltonian of the "fast" modes (see formulae (|TBp and (|TBp). We also 
discuss a possibility to consider quantum chaos as a source of the holographic feature of 
the M(atrix) theory. 

A classical dynamical system is defined by its phase space P, the dynamical fiow 5*^ 
and the invariant measure /i. For the Hamiltonian system one considers the reduced phase 
space obtained by fixing the energy and other (if any) integrals of motion. The measure 
/i is the corresponding restriction of the Liouville measure Yldqidpi. The system is said 
to exhibit a chaotic or stochastic behavior if it is ergodic and moreover it is unstable, i.e. 
it has a positive Lyapunov exponent. 

There are different levels of chaos and they can be specified for instance by /('-property, 
central limit theorem, exponential decay of correlations etc. |[T^-[|1^. Chaos is often 
quantified by computing Lyapunov exponents. If the dynamical system is defined by 
means of the system of differential equations Xi = Fi{x) then the Lyapunov exponent x 
of a solution Xj(t) is given by 

X=limilog^, (2) 
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where 

i 

and ai{t) is the solution of the equation a = F' {x)a. 

The equations of motion for the action (|l]) in the = gauge read 

X, = [[X,,X,],X,]. (3) 

By varying over Aq one also gets the constraint 

[X„X,]=0. (4) 

The Hamiltonian for the bosonic part of the action (P reads 

H = Tr{\pf-\[X,,X,f), (5) 

where Pi is the momentum conjugate to the Xi. 

It has been proved in |^ that the Hamiltonian of the supersymmetric matrix model 
has a continuous spectrum starting at zero. This result can be interpreted as a man- 
ifestation of the instability of the supermembranes against deformations into stringlike 
configurations. The possible instability of membranes is already evident from the clas- 
sical consideration because the potential energy has valleys through which certain mem- 
brane configurations can escape to infinity without increasing the mass. For the bosonic 
membrane this classical instability is cured by quantum mechanics: the spectrum of the 
quantum Hamiltonian is actually discrete [1^, However, the theory still become un- 
stable if one introduces supersymmetry p. These effects can be seen in a toy bosonic 



model 18 with the Hamiltonian 



H = \{pI+pI) + \xIxI (6) 



The quantum mechanical Hamiltonian has a discrete spectrum |T^, |T^, |TU|. However, 
its supersymmetric version has a continuous spectrum ||^. One observes that certain 
properties of the classical bosonic system are more similar to the properties of the quantum 
supersymmetric system rather then to the quantum bosonic system. Supersymmetry saves 
classics. 

The dynamical system with the Hamiltonian (^) is very interesting because it exhibits a 
complicated chaotic behavior of trajectories. Classical and quantum chaos in the system 



obtained as the reduction of the Yang-Mills theory has been considered in [|18|, ^ 
T9|,see also |3^. Chaos in the Einstein- Yang-Mills equations was discussed in p2|, 

The potential in has zero-directions, the valleys, = and X2 = 0. In these 
directions the particle can move without changing energy and the presence of these hy- 
perbolic valleys is an origin of stochastic behavior of the particle. A typical trajectory 
is plotted on Fig.l. We see that a movement of a particle is limited by four hyperbolas. 
The particle being positioned in one of the valleys starts to oscillate between two nearest 
hyperbolas and after a number of oscillations it changes the valley. So, during the time 
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large enough one can see the particle in any of four valleys. Therefore we cannot treat the 
dynamics perturbatively just making a linearisation around some fixed value of xi and X2- 
To illustrate this instability we present on Fig. 2 a dependence of x{t) = j ^^g p(t) / p{0) on 
t for some initial data. We see that for t large enough goes to a fixed value, x ~ 0.85. 

Note, that the addition of fermions leads to an appearing of the interaction with a 
fermionic current. The new system also has a tendency to produce a chaotic behavior (cf. 
??). 

The system (|^) is, in fact, a particular case of the M(atrix) model Let us consider 
the gauge group SU (2) and take the following ansatz: 

Xi = xiai, X2 = X2cr2, X3 = X3CT3, (7) 



X4 = ... = Xg = 



where Xa = Xa{t) are real valued functions of time, a = 1,2,3 and o"q, are the Pauli 
matrices. Then the constraint @ is satisfied and eqs. @ are reduced to 

= -2(^ (8) 

Eqs d^) could be obtained from the Hamiltonian 

H = ^{pI+pI + pI) + {xlxl + xlxl + xlxl). (9) 

The system is obtained from (|^) by seting X3 = and rescaling the coupling constant. 
The system (pi) is an example of the dynamical system which exhibits a chaotic motion. 



It has been analyzed numerically in [20|. On Fig.3 we show (xi, X2) projection of a typical 
trajectory of the particle and on Fig.4 we plot the corresponding xi^)- During the time 
large enough one can see the particle in any of six valleys. By comparing Fig. 2 and Fig.4 
we see that in the three-dimensional case tends to a constant value faster then in the 
two-dimensional case. 

The ansatz can also be tested on su(3) 

X, = XiT\ I = X^ = Q, (10) 

where are the su(3) generators. The Lagrangian equations of motion put the constraint 

on Xi 

xl + xl-xl-xy = 0. (11) 



One of the possible solutions of (|TTj) is X4 = Xg, x^ = x-j. It yields the system with six 
degrees of freedom with the Hamiltonian 



^=\Y.v\^-Xp\^pI)^\pI^\ E x\x]^\[x\^x\) E 

^ 1=1 ^ Z Z l<j<j<3 ^ 1<*<3 
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where pi is momentum conjugated to Xj. The estimation of the Lyapunov exponent by the 
computer simulation exhibits the stochastic behaviour of the system. For the particular 
case X2 = = X5 = Xs = the Hamiltonian has the form 



^ = 2^? + 4P4 + 2 + 4^4^;? + (13) 



A typical trajectory for ( [T^ ) is plotted on Fig. 5 and it can be seen that they are located 
in the valley xi = 0. The corresponding Lyapunov exponent is presented on Fig.6. 

Note that the ansatz (|^) deals only with the matrix diagonal in (space, isotopic) indices 
and all other degrees of freedom are frozen. One can think that stochastic behavior is an 
artifact of this ansatz and instability will be cured by taking into account the fluctuation 
of the frozen degrees of freedom. To analyze this possibility we consider the model (|l|) 
with only two matrices and relaxed constraints. A relaxing of the Gauss law means that 
we consider an interaction of the "electric" field with some current. Denote: Xi = 0i and 
^2 = 02; then the action reads 



S= dtTT[-<Pl + -^l + -[4>^,4>2Y]. (14) 




This action describes the dynamics of the classical vacuum moduli space in the 
SUSY Yang-Mills theory |2|. 



The dynamics of "fast" modes in ( |T4|) can be reduced in special cases to the following 
two-dimensional Hamiltonian systems 



1/9 9n ,99 (m — n)'^ (m + n)" , , 



= i(P? + pl) + Ai?il + ^, (16) 

where m and n are some real constants (see Appendix, where Xi = f, X2 = g)- These 
constants describe a contribution from the "slow" modes to the Hamiltonian of the "fast" 
modes. We see that the "slow" modes are separated but they bring a singular term to the 
Hamiltonian of the "fast" modes. Note that an influence of the extra term 1/2x1 ^^e 



chaotic behavior of a two-dimensional system has been discussed in the recent paper . 
Typical trajectories for the Hamiltonian (|15D are plotted on Fig. 7 and for the Hamiltonian 
(0) on Fig.8. 

In these special cases one can recover the constraints by taking m = 0. In the case 
(p^) one comes back to the ansatz (^ with 0:3 = 0. But in the case ( pTSf ) we get an 
extra repulsive potential. This potential describes a repulsive from walls located along 
two diagonals Xi = ±X2 in the Xi,X2 plane. A movement of a particle is limited by 
equipotential lines 

Xx^X2 + 2j^^2 _ ^ const, 

and the particle oscillates in one of four valleys (it chooses one of them according to the 
initial data). In each of four allowed domains the potential still has a direction of the 
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instability. Therefore, the presence of the constraint does not affect the chaotic character 
of trajectories. 

Different approaches have been developed to address the question what is the way 
in which classical chaos manifests itself in the properties of the corresponding quantum 
system ||T^ , 14], |18]-[2C], |2^-[0. The simplest manifestation of classical chaos for a 



quantum system is the nature of spectral fluctuations of the energy levels. It was suggested 
that for chaotic systems the statistical properties of the spectrum should be that of the 
random matrix theory with the Wigner-Dyson distribution 

P{E\AE) = A| AE|" exp[-B{AE)% (17) 

where a > 0, whereas the quantum version of a classically integrable system is described 
by the Poisson distribution 

P{E\AE) = aexp[-b\AE\] (18) 



The crucial difference between (p!7| ) and (|I8|) is the behavior for AE 0. For the 
Wigner-Dyson distribution ([T7|) one has P{E\AE) — if AE — > 0, i.e. the density of 
the energy levels at the small scale vanishes. This is interpreted as the repulsive feature 
of the energy eigenvalues for quantum chaotic system. One can try to relate this feature 
with the holographic feature of M(atrix) theory. The holographic principle follows from 
the general consideration involving the Bekenstein-'t Hooft bound on the entropy of a 
spatial region. A holographic theory contains only degrees of freedom which carry the 
smallest unit of longitudinal momentum |^ . The transverse density of partons is bounded 
to about one parton per transverse Planck area, in this sense the holographic theory is 
repulsive. The partons form a kind of incompressive fluid. The repulsive feature of the 
holographic theory seems to be related with the repulsive feature of the distribution of the 
energy levels for quantum chaotic system. Quantum chaos is a source of the holographic 
feature of M-theory. A connection between the large limit, random matrix theory and 
entropy of black holes is discussed in [^-@. 

In conclusion, we have shown that M(atrix) theory exhibits the classical chaotic motion 
and we have argued that quantum chaos is a source of the holographic feature. The 
statistical properties of the spectrum and other features of quantum chaos in M-theory 
require a further investigation. 
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Appendix 

Here we shall explore the classical dynamics of the system with the action ([1^) and deduce 
the Hamiltonians ( |15|) and (|T6|) . 

The Lagrangian admits two global continuous symmetries. Since the fields 0i and 02 
are the elements of the su(2) algebra the action is invariant under 

0; = f/+0,f/, 0; = [7+0,^7, (A.l) 
U G SU{2). This symmetry yields the conservation of the "angular momentum" 

Mo= [01, 0l] + [02,02]. (A.2) 

Another symmetry corresponds to U(l) transformations which mix the fields 

(j)^ = cos 6 (pi — sin 6 02, 

02 = sin 6' 01 + cos 6* 02. (A. 3) 

It gives one more first integral 

Ar = Tr(0i02-0i02). (A.4) 



The equations of motion for the Lagrangian (111) read 



01= 2A[02, [01, 02]], 02= 2A[0i, [02, 0i]]. (A.5) 
It is convenient to parametrize 0i and 02 as follows: 

01 (t) = V2U+{t) (^/(t) cose(t) - ^g{t) sin^(t)) U{t), 

Mt) = V2U+it) ^^f{t)sme{t) + ^g{t)cose{t)^ U{t), (A.6) 

where f(t), g(t), 0{t) are real functions and U{t) is a SU(2) group element. 

Let us analyze the consistency of the parametrization (|A.6| ). The variables 0i and 02 



could be treated as vectors in the internal isotopic space. At any time they can be rotated 
to belong to some coordinate plane, say (2,3), by using an U{t) G SU{2): 

01 = (0, ^1, ^2), 02 = (0, V-s, ^-4), (A.7) 

that fixes U{t) up to rotation around the 1-axis. This rotation could be used to impose 
the following constraint: 

+ = 0. (A.8) 

The rotation angle to fulfill ( [A.8|) 



tan 2% 



IS 

2(^lV'2 + V'3V'4 



V^l - ^2 + - 
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Two vectors in the plain obeying ( [A .81 ) could be parametrized by three variables /, g, 6. 
Namely : 

f = ^\/i[+i4, 9 = -^M+i4, tan^^ = -^ 

In this coordinate system the Lagrangian acquires the form 

L = ^{f' + g' + if + g')e' + if + Tr - 2i Tr UU^a^fgO 

+ Tr UU+asUU+crs + Tr UU+a2UU+a2 - Xfg^. (A.9) 

Using "fast-slow" modes terminology we can say that / and g are the "fast" modes 
and 9 and U describe the "slow" modes. 

Note that 6* is a cyclic coordinate and the corresponding first integral is just N = n 
with 

N=if + g')e + 2fgh, (A.IO) 

here we denote: UU~^ = I = ^Cjlj. 

The equations of motion for / and g in this parametrization read: 

/= fO' + fill + ll) + 2geh - 2\fg\ (A.ll) 

9= gQ^ + g{l\ + l^) + 2fQh - 2Xfg, (A.12) 

The remaining three equations of motion are nothing but the angular momentum con- 
servation law Mq = 0. It is convenient to introduce M - the angular momentum in the 
moving frame M = UMqU^ with the covariant conservation law 

M+[M,/] = 0. (A. 13) 

The explicit expression for M reads 

M = lai {if + g')h + 2fge) + ia2{f%) + i(y;{g%). (A.14) 
The Lagrangian can be expressed in terms of M and as follows 

L = \{p + g')- Xfg' + ^i9N- Tr /M). (A. 15) 

There is a wide class of solutions with [M, I] = 0, i.e. where M is conserved. It follows 
from eq.( |A.14|) that in this case only one component of / can be non-zero. 



In the first case: I2 = I3 = eqs. ( |7rT0| ), { \KJ^ ) and {\KJ^ give 



2fgh + if + g')e = n 
if + g')h + 2fge = m. 

The matter of simple calculation to verify that the dynamics of the (/, g) system (eqs. 
( [A?Tl| ), ( |7rT^ )) is governed by the "effective" potential 



; dV dV 

The second case: /i = /a = 0, gives 

Ml = 2fge 

and to keep M parallel to / one has to put = 0. In this case the effective potential is 

^ = A/V + ^, (A.18) 

where m = PI2 is the first integral. 

The third case /i = ^2 = is quite similar to the second one with the exchange 2^3. 
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. 1. Typical trajectory of the two-dimensional system with Hamiltonian (§) 
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Fig. 2. Lyapunov exponent of the two-dimensional system with Hamiltonian (Q) 
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3. Typical trajectory of the three-dimensional system with Hamiltonian (^) 
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Fig. 4. Lyapunov exponent of the three-dimensional system with Hamiltonian (|^) 
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5. Typical trajectory of the two-dimensional system with Hamiltonian ([TSl) . 
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